Compared analysis of proton electromagnetic form factors in space-like and time-like 

regions 



E. Tomasi-Gustafsson 

IRFU,SPhN, 91191 Gif-sur-Yvette and CNRS/IN2P3, IPNO, UMR 8608, 91405 Orsay, France 

New possibilities of high precision measurements of hadron form factors in annihilation and scat- 
tering reactions over an unexplored kinematical region suggest a compared analysis, in view of a 
global description of the nucleon structure. 
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I. INTRODUCTION 

G\ • 

The study of hadron electromagnetic form factors (FFs) is a very active field of high energy physics since a few 
^vq \ decades, as FFs are fundamental quantities which contain information on the internal structure of composite particles. 
, ■ They constitute a privileged playground for the test of theoretical models, and should reflect the transition from the 
non perturbative regime, where effective degrees of freedom describe the nucleon structure, to the asymptotic region 
where QCD applies. The possibility to transfer high momenta, and therefore to access small internal distances, allows 
lO ' to test pQCD predictions, such as quark counting rules and helicity conservation. 

In the space-like region, high precision measurements in an extended kinematical range are an important part of 
the present and future experimental program at Jefferson Laboratory (USA). In the time- like region, a program is 
foreseen by the PANDA collaboration at FAIR (Germany), using high intensity antiproton beams up to 15 GeV 
■+7* ■ kinetic energy. Similar studies are also discussed as part of the experimental program at electron positron colliders, 
in Frascati (Italy), Novosibirsk (Russia), Beijing (China). 

The traditional way to measure electromagnetic hadron FFs is based on elastic electron proton scattering + p <-> 
e~ + p and on the annihilation reactions related by crossing symmetry e + + e~ <-> p + p, assuming that the interaction 
occurs through the exchange of one virtual photon, of mass q 2 . In recent years, very surprising results have been 
obtained in ep elastic scattering, due to the possibility of applying the polarization method 1]: the electric and 
magnetic distributions inside a proton do not have the same dipole dependence, as a function of q 2 0], as it was 
■ previously assumed. 

The understanding and the interpretation of the data at large momenta in the full kinematical region requires to 
T^j" . investigate carefully not only the nucleon structure but also the reaction mechanism. The simple extrapolation of 
models and methods should be taken very carefully. A large debate recently arose, due to inconsistencies among form 
, factors extracted from polarized and unpolarized experiments in space-like region (for a review, see Q). As no bias 
has been found in both types of experiment and as the extraction of form factors follows the same formalism (based 
Q\ [ on one-photon exchange), possible explanations are related to higher order radiative corrections. 

The measured observables are the differential cross section in unpolarized ep scattering and the ratio of the longi- 
tudinal to transverse proton polarization in elastic scattering of longitudinally polarized electron on an unpolarized 
proton target. Radiative corrections are very large for the unpolarized cross section, and are neglected in polarization 
experiments. High order corrections have not yet been applied to the data. The presence of two-photon exchange 
would induce a more complicated spin structure of the matrix element and it has been discussed in the frame of a 
compared analysis of space-like 0, Q and time- like data [E 0, HI • 
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II. FORMALISM 



A. Space -like region 



Assuming one-photon exchange the reduced elastic differential cross section for ep elastic scattering, may be written 



as 



m 



ared{e,Q 2 ) = e(l + r) 
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1 + 2— sin 2 (9/2) 



AE 2 sin 4 (9/2) da 
a 2 cos 2 (9/2) dh 



= rG 2 M (Q 2 )+eG 2 E (Q 2 ), 



e = [l + 2(l + T)tan 2 (6>/2)]- 1 , r = Q 2 /(4m 2 ), Q 2 = -t 
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where a = 1/137, m is the proton mass, E and 9 are the incident electron energy and the scattering angle of the 
outgoing electron in the laboratory system, respectively. Gm(Q 2 ) and Ge(Q 2 ) are the magnetic and the electric 
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FIG. 1: Contribution of the Ge(Q 2 ) dependent term to the reduced cross section (in percent) for e = 0.2 (solid line), e = 0.5 
(dashed line), e = 0.8 (dash-dotted line), in the hypothesis of FF scaling (thin lines) or following Eq. (|3} (thick lines). 



proton FFs, functions of Q 2 , only. Measurements of the elastic differential cross section at different angles for a fixed 
value of Q 2 allow Ge(Q 2 ) and Gm{Q 2 ) to be determined as the slope and the intercept, respectively, from the linear e 
dependence @ . The normalization is chosen in order to have static values proportional to the proton electric charge 
and magnetic moment fi, respectively Ge(0) = 1 and Gm(0) = fi. 

From unpolarized cross section measurements the determination of Ge and Gm has been done up to Q 2 ~ 8.8 GeV 2 
[TTI | and Gm(Q 2 ) has been extracted up to Q 2 ~ 31 GeV 2 [l2j under the assumption that Ge — 0, and it is often 
approximated, for practical purposes, according to a dipole form: Gd(Q 2 ) = (1 + Q 2 /0.71 GeV 2 ) -2 . Polarization 
transfer measurements suggest a monotonical decrease of the ratio /j,Ge(Q 2 )/Gm(Q 2 ) with Q 2 : 

[iGeIGm = 1 for Q 2 < 0.4 [(GeV/c) 2 ], ^G E /G M = 1.0587- 0.14265 for Q 2 < 6 [(GeV/c) 2 ], (3) 

at larger Q 2 , at least up to Q 2 ~ 6 GeV 2 , deviating from unity as Q 2 increases Q. 

At large Q 2 the contribution of the electric term to the cross section becomes very small, as the magnetic part is 
amplified by the kinematical factor r. This is illustrated in Fig. [1] where the ratio of the electric part to the reduced 
cross section Fe = eG 2 E (Q 2 ) /a re d(0, Q 2 ), is shown as a function of Q 2 , for different values of e. The different curves 
correspond to different values of e, assuming FFs scaling (thin lines) or in the hypothesis of the linear dependence of 
Eq. ([3]) (thick lines). In the second case, one can see that, for example, for e = 0.2 the electric contribution becomes 
lower than 3% starting from 2 GeV 2 . This number should be compared with the absolute uncertainty of the cross 
section measurement. When this contribution is larger or is of the same order, the sensitivity of the measurement to 
the electric term is lost and the extraction of Ge(Q 2 ) becomes meaningless. 

Higher precision can be obtained in polarization experiments. As it was firstly shown in [l[, measuring the polar- 
ization of the outgoing proton in the scattering of longitudinally polarized electrons on an unpolarized proton target, 
gives access to an interference term which contains the product Ge(Q 2 )Gm{Q 2 ) and it is more sensitive to a small 
electric contribution than the cross section itself. The following expressions hold for the transverse and longitudinal 
components Pt and Pg of the proton polarization vector (in the scattering plane) in terms of the proton electromagnetic 
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FFs: 

DP t = -2\cotL r^G E G M ,DP t = A^±^,/^G 2 M , (4) 
2 V 1 + r m \ 1 + t 

£> = 2rG M + cot — . (5) 

Z 1 + T 

where E' is the scattered electron energy and D is proportional to the differential cross section with unpolarized 
particles. So, for the ratio of these components one can find the following formula: 

Pt _ . 9 m G E (q 2 ) 

pr- 2cot 2ET^G^m (6) 

which shows the direct link between the polarization components of the recoil proton and the electric and magnetic 
FFs. 

The results obtained with such technique display a large precision compared to the Rosenbluth separation, due to 
the large sensitivity to the electric FF. Moreover, the electron beam polarization as well as the analyzing powers of 
the polarimeter cancel in the ratio, reducing the systematic errors. 



B. Time-like region 

Due to unitarity, in TL region hadron FFs are complex functions of q 2 . The unpolarized cross section depends on 
their moduli, the measurement of which is, in principle, simpler than in SL region where the Rosenbluth separation 
requires at least two measurements at fixed q 2 and different angles implying a change of incident energy and scattered 
electron angle at each q 2 point. In TL region, the individual determination of \Ge \ and \Gm \ requires the measurement 
of the angular distribution of the outgoing leptons, at fixed total energy s = q 2 . All the information, of the nucleon 
structure as well as of the reaction mechanism, as discussed below, is contained in the differential cross section. 

The differential cross section for the annihilation process 

P + P ^£ + + r, £^e,n (7) 

was first obtained in ref. [HI]. It is expressed in terms of the proton electromagnetic FFs as: 

da ira 2 



d(cOS0) %m 2 T^T{T - 1) 



[T|G M | 2 (l + cos 2 0) + |G E | 2 sin 2 #] , (8) 



where 9 is the electron production angle in the center of mass system (CM). The cos 2 9 dependence of Eq. (JSj) results 
directly from the one-photon exchange mechanism, assuming that the spin of the photon is equal to one and that 
the electromagnetic hadron interaction satisfies C-invariance. This corresponds, by crossing symmetry, to the linear 
Rosenbluth cot 2 (9/2) dependence fl4| . 

The electric term is accompanied by a dependence in sin 2 9. It means that, whatever is the model used for |Gb(9 2 )| 2 , 
it has a maximum at cos 9 = and vanishes at con 9 = ±1. The magnetic term has a maximum at cos 9 = ±1, which 
equals to 2t|Ga/(<? 2 )| 2 and a minimum at cos 9 = 0, which equals to t|Gm(? 2 )| 2 - The relative contribution of the 
electric a E (dashed lines) and magnetic <jm (solid lines) terms to the differential cross section are illustrated in Fig. 
[5k., for two values of q 2 = 5 and 8 GeV 2 . 

The total cross section is a quadratic combination of FFs, which does not contain interference terms: 

(2t\G m \ 2 + \G e \ 2 ). (9) 



6m 2 T-\/ t(t — 1) 

Let us introduce the angular asymmetry, A, which enhances the different angular behavior of the electric and magnetic 
terms with respect to 9 = 90°. One can express the angular dependence of the differential cross section in terms of 
the angular asymmetry A as: 

-^- = ( t [1 + .4cos 2 #], (10) 
a(cos 9) L J 

where ctq, the differential cross section at 9 = tt/2, and A can be written as functions of the FFs as: 



1,„ ,o\ . t\Gm\ 2 -\G r \ 2 r-TZ 2 „ \G f 



a ° = vv~ (J Gm|2 + t IGe1 ) 5 A = t\gL? + Kl 2 = TTt^' n = wk- 
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FIG. 2: (a) Relative contribution of the electric <je (dashed lines) and magnetic om (solid lines) terms to the differential cross 
section for p + p — > e + + e~ , as functions of cos# for two different values of q 2 : 5 GeV 2 (red lines) and q 2 =8 GeV 2 (black 
lines) ; (b) Relative contribution of the electric (dashed lines) and magnetic (solid lines) terms to the total cross section Fe and 
Fm (black lines) and to the angular asymmetry, Ae and Am (red lines). 



The angular asymmetry A lies in the range — 1 < A < 1. For Ge = one obtains A = 1 and for Ge — Gm one 
obtains A = (r - l)/(r + 1). 

The electric and magnetic contributions to the total cross section and to the angular asymmetry are illustrated in 
Fig. [2Jd, as function of q 2 . The unphysical region is indicated by a dashed area. Although the magnetic contribution 
largely dominates in the physical region, the relative contribution of the electric term to the cross section is larger 
than 10% for q 2 < 15 GeV 2 , and it is even larger for the angular asymmetry. 

From the total cross section, it is possible to extract \Gm \ under a definite hypothesis on the ratio. The experimental 
results are usually given in terms of |Gm|, under the hypothesis that Ge = or \Ge\ = \Gm\- The first hypothesis 
is arbitrary whereas the second one is strictly valid at threshold only, and there is no theoretical argument which 
justifies its validity at any other momentum transfer, where q 2 ^ Am 2 . However, Ge plays a minor role in the cross 
section and different hypothesis for \Ge\ do not affect strongly the extracted values of Gm, due to the kinematical 
factor r, which weights the magnetic contribution to the differential cross section and makes the contribution of the 
electric FF to the cross section smaller and smaller as q 2 increases. 

The individual determination of the FFs in time-like region has not yet been done. The ratio R = Ge/Gm has been 
determined from a two parameter fit of the differential cross section, by PS170 at LEAR [15j, and more recently by 
the BABAR Collaboration using initial state radiation, e + + e~ -^p + p + -f 16]. Data are very limited and affected 
by large errors, mainly due to statistics. The results from BABAR suggest a ratio larger than one, in a wide region 
above threshold, whereas data from [15] suggest smaller values. The results from PANDA are expected to clarify this 
issue. 



C. Discussion 

A general illustration of the world data for the proton form factors is given in Fig. [3l where proton FFs are shown 
as function of \q |, allowing a straightforward comparison in the whole kinematical region. In order to eliminate the 
steep q 4 dependence, all FFs are rescaled by the dipole function. 

From top to bottom, one can see the magnetic proton FF in TL region, under the assumption \Ge\ = |Gm|, the 
magnetic proton FF in SL region which is obtained for q 2 < 8.8 GeV 2 obtained under the assumption \Ge \ = (blue 
circles) and the electric FF in SL region. Two series of data clearly show the discrepancy between unpolarized (red 
triangles) and polarized (green stars) measurements. 

The expected precision of the future measurements with PANDA (black solid squares) is shown in comparison with 
the existing data. For PANDA each point corresponds to an integrated luminosity of 2 fb _1 , which can be obtained 
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FIG. 3: World data on proton form factors, in time and space-like regions, as functions of |g 2 |, rescaled by dipole. From top to 
bottom, magnetic FF in time-like region including PANDA simulated results (black solid squares), magnetic FF in space-like 
region (blue circles), electric FF in space-like region, from unpolarized (red triangles) and polarized (green stars) experiments. 



in four months of data taking. These results have been obtained in frame of Montecarlo simulations, which takes into 
account the geometry of the detector, efficiency and acceptance and is based on a realistic parametrization of FFs 
One can see that PANDA will cover a large kinematical range and bring useful information with respect to the 
determination of the asymptotic region. 

A more direct representation of FFs is given by the Dirac F\ and Pauli Fi FFs, which are linear combinations of 
Ge and Gm- PQCD predicts the asymptotic behavior F\ ~ Q~ 4 and F^ ~ Q~ 6 which is followed by the Rosenbluth 
measurements, but not compatible with polarization data, which suggest instead the following ratio: F2/F1 ~ 

The values of Gm in the TL region, obtained under the assumption \Ge\ = \Gm\, are larger than the corresponding 
SL values. A difference up to a factor of two in the absolute values in SL and TL regions can be seen also for other 
hadron FFs, including pions and neutrons, up to the largest value at which TL FFs have been measured. This has 
been considered as a proof of the non applicability of the Phragmen-Lindelof theorem, or as an evidence that the 
asymptotic regime is not reached [l||. The Phragmen-Lindelof theorem constrains definitely FFs in TL and in SL 
regions to have the same value at large q 2 . This theorem has other applications in particle physics, such as, for 
example, the Pomeranchuk theorem, concerning the asymptotic behavior of the total cross sections for a + b and a + b 
collisions (a and b any hadrons): aT(ab) = ax{ob). However, to be rigorous, the applicability of this theorem to FFs, 
which seems evident, has not yet been proved. 

In principle, asymptotic properties should be discussed for F\ and F%. The analyticity of FFs allows to apply the 
Phragmen-Lindelof theorem which gives a rigorous prescription for the asymptotic behavior of analytical functions: 

lim F^ SL \q 2 ) = lim F^ TL \q 2 ). (12) 

q 2 — > — 00 q 2 — >oo 

This means that, asymptotically, FFs have the following constraints: 1) the imaginary part of FFs, in TL region, 
vanishes: ImFi(q 2 ) — > 0, as q 2 — > 00; 2) the real part of FFs, in TL region, coincides with the corresponding value in 
SL region, because FFs are real functions in SL region, due to the hermiticity of the corresponding electromagnetic 
Hamiltonian. 

Unfortunately, this theorem does not allow to indicate the physical value of q 2 , starting from which it is working 
at some level of precision. For this aim one needs some additional dynamical information. The assumption of the 
analyticity of FFs allows to connect the nucleon FFs in SL and in TL regions and to extend a parametrization of FFs 
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available in one kinematical region to the other kinematical region. Dispersion relation approaches, which are based 
essentially on the analytical properties of nucleon electromagnetic FFs, can be considered a powerful tool for the 
description of the q 2 behavior of FFs in the entire kinematical region. The vector meson dominance (VDM) models, 
can be also extrapolated from the SL region to the TL region (see Q and Refs. therein). The quark-gluon string 
model [18J allowed firstly to find the q 2 dependence of the electromagnetic FFs in TL region, in a definite analytical 
form, which can be continued in the SL region. 

In order to test these requirements, the knowledge of the differential cross section for e+ + e~ <-> p + p is not 
sufficient, and polarization phenomena have to be studied. In this respect, T-odd polarization observables, which 
are determined by ImF^F^, are especially interesting. The simplest of these observables is the P y component of the 
proton polarization in e + + e~ — ► p + p that in general does not vanish, even in collisions of unpolarized leptons, or 
the asymmetry of leptons produced in p + p — > e + + e~, in the collision of unpolarized antiprotons with a transversally 
polarized proton target (or in the collision of transversally polarized antiprotons on an unpolarized proton target) 
[1> [HI- These observables are especially sensitive to the different parametrizations of FFs, and suggest that the 
corresponding asymptotics are very far [201 ] . 



III. TWO PHOTON EXCHANGE 



As stressed in the introduction, the expressions of the cross section Eqs. ([21 [SJ assume one-photon exchange. In 
principle, the interaction can occur through two (or more) photon exchange (TPE). Although if such mechanisms 
are suppressed by powers of a they could play a role at large q 2 , due to possible enhancement from a mechanism 
where the momentum is equally shared between the photons. In such case the decrease of the cross section due to 
a counting would be compensated by the steep decrease of FFs with q 2 . Recently, the possibility of a sizable TPE 
contribution has been discussed as possible solution to discrepancies between experimental data, on elastic electron 
deuteron scattering [l4[ and elastic electron proton scattering [3[. 

The model independent analysis of experimental observables taking into account the TPE contribution, for ep 
scattering and for the crossed annihilation channels can be found in Ref. [!, H, 0] • The presence of TPE induces a more 
complicated spin structure of the matrix amplitude. In the scattering channel, instead that two real FFs, functions 
of one kinematical variable, q 2 , one has to determine three FFs, complex functions of two kinematical variables, and 
the e linearity of the Rosenbluth formula does not hold anymore. However, it is still possible to measure the real FFs, 
using electron and positron scattering on proton, in the same kinematical conditions, or measuring three T-odd or five 
T-even polarization observables. In the annihilation channel TPE induces four new terms in the angular distribution, 
of the order of a compared to the dominant contribution and which are odd in cos 9. 

Therefore, the non linearity of the Rosenbluth fit in the scattering channel d the presence of odd cos 8 terms in the 
annihilation channel can be considered as a model independent signatures of TPE (more exactly, of the real part of 
the interference between one and two photon exchange). Evidence of TPE, based on these signatures has not been 
found in the experimental data on electron elastic scattering on particles with spin zero [2l| , one half [ij , and one [3| • 
An analysis of the B ABAR data [l(| also does not show the evidence of two photon contribution 0] . On the basis of 
simulation studies, it can be shown that the future PANDA experiment will be sensitive to a TPE contribution > 5% 
of the main (one photon) contribution [22| . 

Let us stress that the main advantage of the search of TPE in TL region is that the information is fully contained 
in the angular distribution. TPE effects cancel (are singled out) in the sum (difference) of the cross section at 
complementary angles, allowing to extract the moduli of the true FFs [1,0]. TPE effects also cancel if one does not 
measure the charge of the outgoing lepton. 



IV. CONCLUSION 



In the next future the knowledge of electromagnetic proton FFs will be extended in a wide kinematical region, 
allowing a unified description in both SL and TL regions. It will be possible to clarify both issues, the reaction 
mechanism and the proton electromagnetic structure at short distances. In particular, the individual determination 
of FFs will be possible in TL region, for moderate q 2 values. These data are expected to constrain nucleon models. 
For q 2 > 20 GeV 2 , where the electric contribution should become negligible, the validity of asymptotic properties 
predicted by QCD will be tested. Note that there is no principal reason for which the electric and the magnetic FFs 
should reach the same asymptotics, at the same q 2 . 

Due to crossing symmetry properties, the reaction mechanism should be the same in SL and TL regions, at similar 
values of the transferred momentum. If TPE is the reason of the discrepancy between the polarized and unpolarized 
FFs measurements in SL region, a contribution of 5% is necessary to bring the data in agreement in the \q 2 \ range 
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between 1 and 6 (GeV/c) 2 . Such level of contribution will be detectable in the PANDA experiment. In Ref. [23[ the 
discrepancy has been attributed to the method of calculating radiative corrections. Radiative corrections are specific 
to each of these reactions, therefore a comparison of the data issued from the three channels, ep scattering, e + e~ and 
pp annihilation, will shed light on the reaction mechanism. 

This work has been initiated within a long term collaboration with Prof. M. P. Rekalo and it is largely based 
on common work with Dr. G.I. Gakh. The members of the PANDA Group at IPN Orsay are acknowledged for 
interesting discussions and remarks. 
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